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D . Abstract 



p. 
D 

We present a new analysis of the rare decay Kl —>■ vr^e+e" taking into account im- 
^ ■ portant experimental progress that has recently been achieved in measuring Kl — * n^'yy 
. and Ks — > TT^e~^e~. This includes a brief review of the direct CP- violating component, 
a calculation of the indirect CP-violating contribution, which is now possible after the 
measurement of Ks 7r^e~^e~, and a re-analysis of the CP conserving part. The lat- 
ter is shown to be negligible, based on experimental input from Kl 7r°77, a more 
general treatment of the form factor entering the dispersive contribution, and on a com- 
parison with the CP violating rate, which can now be estimated reliably. We predict 
B{Kl 7r°e+e-) = (3.2tJ:|) x 10"^^ in the Standard Model, dominated by CP violation 
with a sizable contribution (~ 40%) from the direct effect, largely through interference 
with the indirect one. Methods to deal with the severe backgrounds for Kl —>■ 7r°e+e~ 
using Dalitz-plot analysis and time- dependent K^-Ks interference are also briefly dis- 
cussed. 



1 Introduction 



The flavour- changing neutral-current transition Kl ir e~^e~ has been recognized since 
a long time to be one of the most interesting rare kaon decays. It shows an intriguing 
interplay of short and long distances, leading to a sum of comparable CP-conserving, 
direct- and indirect-CP-violating contributions [r. Until very recently, it was impossible 
to estimate all these contributions with good accuracy, or to predict the total Kl — > 
7r°e^e~ rate. As a consequence, it was not clear to which extent this mode could be 
used as a probe of the Cabibbo-Kobayashi-Maskawa mechanism of CP violation. The 
situation has now been changed substantially by two new experimental results of the 
NA48 Collaboration: the observation of the Ks —>■ vr^e^e" decay j2] and the precise 
measurement of the —>■ n^'j'j spectrum at small diphoton-invariant mass [3]. The 
former allows us to evaluate the indirect-CP-violating part of the amplitude, whereas the 
latter help us to estimate the CP-conserving contribution. The purpose of this paper 
is a complete reanalysis of the Kl 7r^e^e~ decay in view of this new experimental 
information. 

The use of experimental data to estimate indirect-CP-violating and CP-conserving 
contributions of — > n^e'^e' is not completely straightforward. The most delicate 
issue concerns the CP-conserving contribution and, particularly, the dispersive part of 
the Kl — > n^Yl* ~^ 7r°e+e~ amplitude. So far, this part of the amplitude has been 
estimated employing a specific model-dependent ansatz for the behavior of the Kl 
7]-o^*^* vertex with off-shell photons |3]. Here we generalize previous analyses using a 
more general parameterization of the Kl — > tc^'J*'-/* form factor, which satisfies both low- 
and high-energy constraints and helps us to estimate the theoretical uncertainties in this 
calculation. Moreover, we show how to extract the information on the on-shell Kl — > 
7r°77 amplitude relevant to Kl — > n^e~^e~ in a model-independent way, without relying 
on mo del- dependent assumptions on the former, such as the vector-meson-dominance 
parameterization in terms of ay- As a result, we are able to derive a conservative upper 
bound on the total CP-conserving contribution of the Kl — > n^e^e' rate, which turns 
out to be well below the level of the interesting short- distance component. 

As far as the CP- violating amplitude is concerned, the NA48 result on B{Ks 
7r°e+e~) provides us with an unambiguous indication that the indirect-CP- violating con- 
tribution is large and cannot be neglected. Here the most delicate issue is the model- 
dependent sign of the interference between direct- and indirect-CP-violating components 
of the amplitude. As we shall show, the measured value of B{Ks — > 7r°e"'"e~), together 
with theoretical arguments both of perturbative and non-perturbative nature, provides a 
good indication in favour of a positive interference. Following this indication, we predict 
B{Kl — > vr°e"''e~)sM ~ 3 x 10^^^, with a negligible CP-conserving component, about 
40% due to the clean short- distance direct-CP-violating amplitude (mainly through the 
interference with the indirect-CP-violating one) and an overall theoretical error that to a 
large extent scales with the experimental error on B{Ks — > 7r°e"'"e~). 

From a purely theoretical perspective we thus conclude that Kl —>■ vr^e^e" appears 
as one of the most interesting candidates for precision tests of CP violation in AS = 1 
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transitions. The most serious problem to reach this goal is the so-called Greenlee back- 
ground ISlini, that is, the large irreducible experimental background induced by the decay 
Kl — ^ 77e"'"e~. The enhancement of the total Kl n^e^e^ rate due to the positive 
interference between direct- and indirect- CP- violating amplitudes provides good news in 
this respect, suggesting that future high-statistics experiments could be able to detect 
the Kl 7r°e"'"e~ signal over Greenlee's background. As we shall show, a Dalitz-Plot 
analysis and, especially, time-dependent measurements could provide additional handles 
against this problem. 

The paper is organized as follows: in Section El we briefly review the prediction for the 
short-distance direct-CP-violating amplitude of — > TT^e~^e~ within the Standard Model 
(SM). Section 121 contains one of the two main results of this work, namely the new conser- 
vative estimate of the CP-conserving branching ratio. The second main result is presented 
in Section m where we analyse the interference between direct- and indirect-CP- violating 
amplitudes, and the prediction for the total rate. A discussion of possible Dalitz-Plot 
and time-dependent analyses against the CP-conserving and, especially, Greenlee's back- 
ground is presented in Sections El and (HI The results are summarized in the Conclusions. 
The Appendix contains a self-contained model-independent analysis of the Kl 7r°77 
amplitude in chiral perturbation theory (CHPT) beyond the lowest non-trivial order. 



2 Short-distance contribution to Kl n^e^e 

The direct-CP- violating transition K2 TT^{e~^e~)j=i, where the lepton pair forms a 
vector or axial-vector state, is dominated by short- distance dynamics and is calculable 
with high accuracy in perturbation theory PflEj. Following the notation of 8j, the effective 
Hamiltonian necessary to compute this amplitude at next-to-leading order accuracy can 
be written as 

^ 6,7V 

n[)f=' = V:^V^,[ + Ty,ifi))Q,ifi) + Ty,A{Mw)Q7A{Mw) ] + h.c. (1) 

V ^ i=l 

where r = —(y^gVtd)/{y*gVud) and Vij denote CKM matrix elements. Here Q12 are the 
current-current operators, Q3...6 the QCD penguin operators, and 

Qjv = ^7^(1 - 75)^^7/ , Q7A = sr{l - l^)dh,id (2) 

the leading electroweak operators. Employing the standard CKM phase convention, the 
overall factor V*gVud is real and direct CP violation is induced only by terms proportional 
to Q'r. Since yi = 1/2 = 0, only Qyy, Qja and, in principle, the QCD penguin operators 
Q3...6, are relevant to estimate the K2 7r°(e+e^)j=i amplitude. 

Ignoring for the moment the contribution of Q3...6, whose matrix elements vanish at 
the tree-level, one has 

A{KL^ne+e-)^^ = - {^{p^)e+ {k,)e- {k2)\n[)f='\K2{p)) 
G 

= i^^{V*Vtd)f+iz){p^+p)''[y7vu{k2h^,v{ki) +y7Auik2hi^75vih)] ■ (3) 
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where f+{z = /m\) f» 1 + z{mK /frtpY is the form factor of the vector current and 
= {p — Pn)^- If considered alone, this theoretically clean part of the amplitude leads to 



B{Kl IT e+e )cpv-dir = ^(;^+) |2 (^7A + yiv) mVtsVtd) 



(2.4 ±0.2) X 10 



-12 



10- 



(4) 



where 



= '^{V*Vtd) ^ (1.36 ± 0.12) x 10"^ (5) 

and the numerical coefficient in (jH) has been obtained using Vus = 0.2240 ± 0.0036 0, 
a{Mz) = 1/129, yjA = -(0.68±0.03)a(Mz) andyjv = (0.73±0.04)a(Mz), corresponding 
to mt{mt) = 167 ± 5 GeV and a renormalization scale for |/7y chosen between 0.8 and 
1.2 GeV IE]. 

The contribution of the QCD penguin operators cannot be estimated with good accu- 
racy, due to unknown hadronic matrix elements. However, Q3„s are expected to yield a 
negligible correction to the leading electroweak operators: 

^y,(/i)(7re+e-|Q,|i^) « yrv{7ie+e-\Qjv\K) . (6) 

i=3 

This assumption is strongly supported by i) the corresponding relation for the quark-level 
matrix elements 

6 

^?/i(/i)(c/e+e"|(5i|s) < ?/7v((ie+e"|Q7v|s) , (7) 

which can be checked explicitly in perturbation theory; ii) the following order-of-magni- 
tude estimate for the impact of the Qi at the hadronic level: 

6B{Kl - 7r°e+e-)Q, ~ ^^^B{K^ vr+e+e") < 10"^^ . (8) 

Given these considerations, Eqs. and (jH) can be considered as precise estimates of 
the Kl — >■ 7r^e~^e~ direct-CP-violating amplitude, and its corresponding branching ratio, 
within the SM. 

It is worth to recall that, being dominated by short distances, this part of the Kl — >■ 
n^e~^e~ amplitude is strongly sensitive to possible non-standard contributions. In partic- 
ular, axial and vector components of the amplitude could separately or both be enhanced 
by factors of C(l) with respect to the SM case. In optimistic but still realistic super- 
symmetric scenarios, the short- distance branching ratio in Eq. could exceed the 10~^^ 
level [ini. 
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3 Estimates of the CP-conserving amplitude 



The \7T^e~^e~) state is not, in general, a CP eigenstate: its transformation under CP 
depends on the angular momentum of the lepton pair. The short-distance Hamiltonian in 
Eq. (P) creates the lepton pair in a state of total angular momentum J = 1 (and orbital 
angular momentum L = 0) so that |vr''e"'"e~)s.d. has opposite CP with respect to \K2). On 
the contrary, the long-distance process K2 7r°77 —* 7r°e~*"e~ can lead to final states with 
even J, allowed in the limit of exact CP symmetry El IS El El Ej • The two- 
photon exchange is not the only source of CP-conserving contributions to Ki vr^e^e", 
in principle also dimension-8 operators generated by W-hox diagrams can lead to a final 
state with J 7^ 1; however, such contributions turn out to be completely neghgible [T7] . 

In the limit of exact CP symmetry, the most general amplitude describing the [p) 
'^^{P3)i{Qij ^1)1 {(12,^2) transition, with on-shell or off-shell photons, can be written as 



^lfM^2u 



H o {p-qiq2P +p-q2P^qi-p-qip-q2 -qi-q2P^Pj 



m 



K 



where 



y 



p-{qi -q2) 



{qi + q2f 



IGsl = 9.0 X lO"*' GeV-2 . 



m 



K 



m 



(9) 



(10) 



K 



and Gg = GF\VusVud\g?,/ with l^fgl = 5.1. Due to Bose symmetry the invariant ampli- 
tudes A and B are symmetric under the exchange qi ^ q2- 

Restricting the attention to the on-shell case {qf = Q'l = 0), the allowed range of the 
dimensionless variables is 



0<\y\<-\'^'il,rl,z) 



where 



X{a,b,c) =a^ + b^ + c^ -2{ 
and the unpolarized differential rate reads 



dydz 



z^\A + B\ 



0<z<(l-rA 



ac + be) 



rriK 



y 



4 



A(l,r^,^) \B 



(11) 
(12) 

(13) 



In order to obtain the total rate from (fT^ . the integration is to be performed over positive 
values of y only. 

Instead of using the A-B basis, one can employ the S-B amplitude basis, where 
S = A + B. As can be seen from Eq. (IT!^ . in the S-B basis the two amplitudes do not 
interfere in the on-shell unpolarized rate. This indicates that S and B describe transitions 
to different final states. Indeed, neglecting the y dependence, S describes the decay into 
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|77)j=o and B the decay into |77)j=2. Since the transition [77)7=0 k^e') is forbidden 
in the hmit rrie — 0, if S" does not depend on y, only the B amphtude can induce a 
non-negUgible CP-conserving contribution to Kl n^e^e' IT^ IT^ . The possible y 
dependence of S is strongly suppressed by the combined action of Bose symmetry, which 
forbids linear terms, and chiral symmetry. For this reason we shall concentrate in the 
following mainly on the CP-conserving contribution to 7r^e~^e~ generated by the B 

amplitude. 

Within CHPT the first non-vanishing contribution to the on-shell — > 77^77 transi- 
tion is generated at At this order only a y-independent S amplitude turns out to 
be different from zero ^Hj- However, this parameter-free prediction underestimates the 
observed Kl — >■ 71^77 rate roughly by a factor of 2, indicating the presence of sizable 0{p^) 
contributions. These have been widely discussed in the literature jTH [13 CHI UH 1201 121] 
and include both local counterterms (presumably generated by resonance exchange) and 
non-local terms due to unitarity corrections. By power counting, 0{p^) local countert- 
erms do not induce y-dependent terms in S and B. Unitarity corrections do induce a 
dependence on y"^ in S; however, this is numerically rather suppressed 

An updated discussion about the structure of the Kl — > 77^77 amplitude beyond 
C(p^) can be found in the Appendix. As far as the B amplitude is concerned, here we 
simply note that B{z) is well approximated by -B(O), whose value depended on a specific 
combination of 0{p^) counterterms. The most efficient and model-independent strategy 
to determine the magnitude of |-B(0)| is provided by a measurement of B{Kl 7r°77) in 
the low diphoton invariant mass region. Neglecting the kinematical dependence of S and 
B (expected to be very mild in both cases below the vrvr threshold) and setting the cut 
M^^ < 110 MeV, to make contact with the NA48 analysis we find 



B(Kl ^ 7r%7)A/,,< 



110 MeV 



2.1 X lO''' X / dz 
Jz<om9 

2.0 X 10'^ X |5fO)P X 



.^A(l,r^,z)^/^|5(0)r + lA(l,r^,zp|i?(0)p 
l + 0.04.l^(°)l^ 



mow 



(14) 



If S{0)/B{0) = C(l), the experimental ratio 



D _ B{Kl 7r°77)M-,^<110 MeV .^^s 

^~ 2.0 X 10-9 ^ ^ 

provides an excellent approximation to |i?(0)p; moreover, the inequality |S(0)p < Rb 
holds independently from any assumption on S{0)/ B{0). As shown in Fig. Q the ap- 
proximation B{z) ^ B{0) works very well for |-B(0)| > 1, which is the case of interest 
for Kl — > 7r°e+e'. Thus Rb provides a powerful model-independent tool to estimate (or 
constrain) the size of |i?(-2)| in the whole physical range. 

Having determined the absolute value of the B amplitude at qf = q2 = from ex- 
periments, we can use it as a (real) effective coupling for the Kl ^^^l {11)1 {12) vertex 
in Fig. 121 in order to estimate the CP-conserving Kl — ^ n^e^e' amplitude. If we ignore 
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B{z)\ 
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Figure 1: Absolute value of the B amplitude for different values of the O(p^) local coun- 
terterms (see the Appendix). 

the dependence on q\ ^ the two-photon dispersive integral turns out to be logarithmically 
divergent. To regularize it we employ the following ansatz 

B{z,y-ql^^) = B{z)^j{ql<g) (16) 

where the form factor 

4) = l + a + ^i^) + b-, , (17) 

is defined in analogy with the analysis of the Kl — 77 ^ amphtude presented in 

[22] • This structure is dictated by the assumption that VMD plays a crucial role in the 
matching between short and long distances (in the numerical analysis my is identified 
with mp ~ 770 MeV). In order to obtain an ultraviolet convergent integral, as it is the 
case within the full theory, we need to impose the condition 

1 + 2a + 6 = . (18) 

In this way we are left with a single free parameter, that we can vary to estimate the 
theoretical uncertainties of this approach.^ Obviously Eq. ()16|) does not represent the 

^ In principle it could be possible to obtain additional constraints on this form-factor by means of 
Kl — > 7r'^Z+^~7 data |23| : however, the present experimental information is not accurate enough to extract 
any significant constraint j24| . 
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Figure 2: Two-photon contribution to the 



amphtude. 



most general parameterization of the B amphtude off shell. However, we are not interested 
in the detailed structure of this amplitude, rather in its weighted integral relevant to 
Kl —>■ TT^e~^e~ . In this respect the possibility to vary one parameter, in order to test the 
stability of the result, represents an important improvement with respect to the choice 
made in |3]. The latter is recovered as a special case for b = —a = 1. 

Following [11 we can write the CP-conserving Kl{p) — > 7r°(p^)e+(A;i)e^(A;2) amplitude 

as 



M ( ^ 7r° e+ e 



iCPC 



Gsa'^B{z)G{z) 



(19) 



K 



where z = s/mj^ = {ki + k2Y /m\ and G{z) is a dimensionless function resulting from 
the loop integration. Neglecting terms which are suppressed by inverse powers of my and 
eliminating h by means of Eq. (jl8|) . we find: 



G{z) 




1 
9 



4 



'1 



(20) 



The function B{z) is certainly real for z < Am\/m\;. Above the tttt threshold this is no 
longer true; however, as discussed above, if |-B(-2)| ~ 1 the local contribution is dominant 
and, to a good accuracy, B{z) can be approximated by a real constant in the entire phase 
space. For this reason the product |-B(2;)| x '^G{z) (where ^G{z) is obtained from 
by the prescription ln(— m^) = Inm^ -|- m) represents to a good approximation the full 
absorptive contribution of the CP-conserving 7c^e^e~ amplitude. 

Our result for this absorptive part is consistent with the one in jl3lin] but not with the 
result of |4|. We disagree with [4 also in the dispersive part, 1-8(^)1 x ^G{z), computed for 
a = —1, where our form factor is identical to the one of The main difference between 
our result and the one in |lj is that we do not find any singularity in the limit rUe 0. 
The lacking of this kind of singularity was already noticed and discussed in jTSj . Indeed 
it is possible to explicitly check that there are no possible bremsstrahlung contributions 
that could cancel such a singular behaviour. 
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Integrating the amplitude ()19|) over the phase space, at fixed 2;, leads to 



8.1 X 10^1=^ X \B{z)\^\{l,rl,zf/^ 



1 + ( ^mz] 



(21) 



As we shall discuss in more details in Section the kinematical factor A^/^ strongly sup- 
presses the high-2 contributions, strengthening the validity of the approximation B{z) ^ 
B{0). In this limit, integrating over the full phase space we obtain 

B{Kl vr°e+e-)cpc = 7.0 x 10"^^ x |S(0)|^ 

X {1 + [1.4 + 1.4(1 + a) + 0.4(1 + a)2] } . (22) 

The term between square brackets in ()22|) represents the model-dependent contribution 
of the dispersive amphtude: we express it as a function of (1 + a), since (a + 1) ~ 0, 
corresponding to the choice of Ref. j3], denotes the most natural possibility. As can 
be noticed, for reasonable values of the form-factor parameters, namely for both a and 
b = —1 — 2a of 0{1), absorptive and dispersive contributions are of the same order. 
Assuming the term between curly brackets in ()22|) to be smaller than 10, which we consider 
a rather conservative hypothesis, and using the inequality |i?(0)p < Rb we can finally 
write 

B{Kl ^ 7r°e-*-e-)cPc < 3.5 x 10"^ x B{Kl ^ A7)Af,,<iio Mev • (23) 
Using the recent experimental input [3] (see the Appendix) 

5(ir^ ^ A7)Af,,<iio Mcv < 0.9 X 10-« at 90%C.L., (24) 

our present estimate of the maximal CP-conserving contribution to Kl TT^e~^e~ reads 

B{Kl 7r°e+e")cpc < 3 x 10"^^ . (25) 

We stress that Eq. (j25j) is a conservative upper bound. A reasonable estimate of B{Kl 
7r°e+e~)cpc; based on the assumption that real and absorptive contributions are equal, 
would lead to values below 10^^^. 

To conclude, we recall that Eq. has been obtained under the assumption |i?(0)|<;l. 
Thus, according to Eq. (fT^ . it makes sense only if B{Kl T^'^1l)M^^f<iio MeV ~ 2 x 
10~^. If experiments would find that B{Kl —> 7r°77)Af^^<iio MeV is below this figure, 
other sub- leading amplitudes could become relevant. In particular, we have neglected the 
contribution induced by y^-dependent terms in S. The latter leads to a CP-conserving 
amplitude 

M{Kl^ n'^e+e-fc-pT ^ ^Pih- k,){p + p^),u{k,)^^v{k,) , (26) 

which is not only suppressed by the smallness of d'^S/dy'^, but also by extra kinematical 
and numerical factors (the coefficient z/12 arises from the calculation with a point-like 
form factor). These contributions become relevant only if B{Kl 7r°e"'"e~)cpc is in the 
10~^^ range and thus small enough compared to the CP-violating term. 
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4 Indirect-CP-violating amplitude and total rate 



The CP-conserving transitions Ks n^i^i^ and n'^£~^£~ , with the lepton pair in a 

vector state, are dominated by the long-distance process K{p) 7r7 — > Tr{pT,)i~^{ki)i~ {k2) 
j2E]- The decay amphtudes can in general be written as 



A(Ki^ntr) = ^^^W,{z){p + p^ru{k2Mh) , (27) 

where Wi{z) are form factors regular at z = {i = ±, S). The latter can be decomposed as 
a sum of a polynomial piece plus a non-analytic term, W^'^{z), generated by the vrvr loop 
and completely determined in terms of the physical K ^ 37t amplitude ^7\. Expanding 
the polynomial term up to 0{p^) we can write 

Wi{z) = Gprnl (a, + b.z) + W^iz) , (28) 

where the real parameters and hi encode local contributions starting at 0{p^) and 
0{p^), respectively. High-precision data on — > 7r^e"'"e~ by BNL-E865 [2H] have been 
successfully fitted using Eq. fOHj) and lead to 



a+ = -0.587 ±0.010 6+ = -0.655 ± 0.044 . (29) 

Unfortunately, chiral symmetry alone does not help to determine the unknown couplings 
as and hs in terms of a+ and 6+. On the other hand, the non- analytic term Wg'^{z) is 
known to be very small, due to the A J = 3/2 suppression of the CP-conserving Ks 37r 
amplitude. As a consequence, the Ks — > 7r°e"'"e~ rate turns out to be dominated by local 
contributions 

B{Ks ^ 7r°e+e") = [o.Ol - 0.76a5 - 0.216^ + 46.5a| + 12.90565 + 1.446|] x 10"^° 
^ 5.2 X 10'^ X al , (30) 

where the second line follows from the assumption bs/as = mj^/rn^, motivated by VMD. 

The first experimental evidence of the Ks —>■ 7i^e~^e~ transition has been announced 
very recently by the NA48/1 Collaboration at CERN. The observation of 7 events in 
a clean signal region (with 0.15 expected background events) leads to the preliminary 
result 0: 

BiKs ^ vr°e+e-)^,,>i65 Mev = {s.OtH ± O.2) x 10"^ , (31) 

which implies 

\as\ = imt'ofi ■ (32) 

This result is in good agreement with the naive chiral counting expectation as = 0{1) j2Zj 
and, within this framework, it is quite compatible with the ranges expected from large- 
Nc |2n] and resonance-saturation arguments jHOl- Interestingly, the central value of the 
NA48/1 result is also in excellent agreement with a very old prediction by Sehgal, based 
only on VMD and A/ = 1/2 rule 
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Taking into account the interference between direct- and indirect-CP-violating com- 
ponents, the full CP-violating contribution to 7T^e~^e~ can be written as 



B{Kl 7r"e+e-)cpv = 10"'^ x 



2 



C*mix i Cint , + Cdir 



10"^ / \ 10" 



(33) 



where the ± depends on the relative sign between short- and long-distance amplitudes, 
and 

C^mix = 10''\e\'^^B{Ks-^7r'e+e-) = {15.7±0.3)\asW 
= 10'^^(^?^ + ^rV) = 2.4±0.2, 

Cint = 2sin0,yCUX^^^^^=^= (6.2±0.3)|a5| , (34) 

'VjA + Viv 

with 0£ = 43.7°. The numerical expressions for Cmix and Cjnt in terms of \as\ are computed 
assuming bs/as = m|^/m^ (or the same form factor for direct- and indirect-CP- violating 
components) and considering only the quadratic terms in (j30|l : the quoted error reflects 
the impact of the linear terms. On the other hand, the l.h.s equations in ()34j) are valid 
independently of any assumption on Ws{z), whose possible difference from the short- 
distance form factor f+{z) is taken into account by the factor 



[/dzA3/2(l,^,r2)|iy5(^)|2]V2^ [/d^A3/2(l,z,r2)|/+(^)|2ii/2 ' 



where r,r = mT^/iriK and A(a,6, c) = a? + IP' + — 2{ab + ac + be). We also remark 
that in evaluating B{Ks — ^ TT^e~^e~), a = 1/137 has been used for the electromagnetic 
coupling. For a given branching ratio, a different choice of a implies, effectively, a different 
normalization of \as\ without changing the estimate of Cmix- 

Given the recent result in ()32|) . we can conclude that, within the SM, the indirect- 
CP-violating amplitude is the largest component of B{Kl —>■ TT^e^e^) and the full CP- 
violating contribution is well above the CP-conserving one. 



Sign of the interference term 

As can be seen from fl33|) . the relative sign between short- and long-distance contributions 
is a crucial ingredient for estimating the sensitivity of future high-statistics experiments 
to perform precision SM tests in this mode. A prediction of this sign requires a better 
understanding of the dynamical origin of the local couplings and 6j. To this purpose, we 
first note that the experimental determination of the ratio 6+ / a+ does not strictly follow 
naive dimensional analysis, which would predict = 0[m'j^/{A7TF.,^)'^] ~ 0.2. Such 

a large ratio between 0{p'^) and 0{p^) counterterms naturally points to the presence of 
large VMD contributions in these channels. Under the rather general hypothesis that the 
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bi terms in Eq. (PHj) are entirely generated by the expansion of a vector-meson propagator, 
we can re- write the polynomial contribution to Wi{s) as 



zm 



nVMD 



^^VMD I nVMD\ , VMD 



m 



K 



m 



V 



(36) 



where af^^^ denotes 2;-independent non-VMD contributions.^ A comparison with the 
experimental results in Eq. ()29p then leads to 



a 



VMD 



K 



-fo^'^P = -1.6 ±0.1 , 



a 



nVMD 



_exp VMD 



1.0 ±0.1 . 



(37) 



The large value of a"^^^ can be justified, in the charged channel, by the presence of sizable 
pion-loop contributions. On the contrary, one would expect a pure VMD contribution 
in the Ks mode, where the pion term is negligible. We thus expect a^^^ = and 



,VMD 



a^"'^ = as, an assumption that justifies the use of the same form factor for direct- and 
indirect-CP- violating components adopted to derive ((221) • 

To make a theoretical prediction for as and, particularly, to determine its sign, we 
need to take an additional step. The contributions to the amplitude (j2ZI) generated by 
the leading short-distance operator in Eq. ((H), namely Qrv, satisfy the AI = 1/2 isospin 
relation 

If this relation is obeyed by the full VMD amplitude, we should expect 



[a 



VMD\ 

S ){Q7V) 



,VMD 



1.6 ±0.1 . 



(39) 



Given the various assumptions behind this prediction, the agreement with the experi- 
mental result in (jH^ is rather encouraging.^ This brings us to the likely hypothesis that 
the VMD part of the Wi{z) is dominated by the CP-conserving contribution induced by 
Qjv, in the limit of a low matching scale {rriK ^ fJ^ ^ my).^ The form factors computed 
under this assumption from the short- distance Hamiltonian cannot be trusted in detail, 
since zjy{jj) exhibits a strong scale dependence (that should be matched by the matrix- 
elements of four-quark operators) and we would like to extrapolate it beyond the validity 
region of perturbation theory. However, the VMD structure of the form factors ensures 
we don't need to push this scale too low. Moreover, we shall employ this procedure only 
to fix the relative sign between direct- and indirect- CP- violating contributions. Using the 



^ A similar generic decomposition between VMD and non-VMD contributions has been successfully 
tested in the Kl 7r+7r^7 mode [521155] . 



The A/ = 1/2 relation between 0(p ) local counterterms of K'^ 



and Ks 



amplitudes had already been adopted in Ref. and later on also in Ref. PO). However, in these 
works no distinction is made between VMD and non-VMD contributions to a+, which turns out to be a 
fundamental ingredient to obtain a phenomenologically acceptable prediction for as- 

Note that, consistently with this hypothesis, the contributions generated by Qjv satisfy the relation 
(6_i_/a+)^Q.j.^,^ = {i>s/cis){Q7v) ~ which follows from the presence of the vector form factor f+{z) 

in Eq. ©. 
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Q-jv part of the short- distance Hamiltonian to determine the full Kl 
form factor, we find 



e vector 



2 4:TTyrv{^M) 



V2i 



a 



e e 



(40) 



Since the ratio Zjvilj) IVivilA is negative for /i < 1 GeV j?! |H], we conclude that the 
most likely possibility is a positive interference between direct- and indirect- CP- violating 
components. 

The prediction of a positive interference between direct- and indirect-CP-violating 
components of Kl Ti^e^e~ had already been reached in Ref. [Tj. This conclusion is 
reinforced by the observation that the perturbative value of as, computed from {QrvifJ^)), 
is smaller than the experimental result in (jH^ for /i = 1 GeV, but it grows, i.e. it goes in 
the right direction, for smaller values of fi. 

The consistency of this prediction can be further confirmed by the following comple- 
mentary reasoning: if we trust the AI = 1/2 + VMD argument leading to Eq. we 
have fully established the sign of as within chiral perturbation theory. This sign follows 
from the sign of a+ which, in turn, is determined experimentally by the interference be- 
tween local and non-local terms in (j2Hl)- We have thus established the sign of as in terms 
of the sign of Gs, the overall couphng of the (8l, l^j) non-leptonic weak chiral Lagrangian 
(see e.g. Ref. j^). This implies that the Kl vr^e^e^ vector form factor can be written 
as 



a5|e|e''^^sgn(G8 



. 47r|/7y(/i) 

I 7= 



V2a 



(41) 



with a positive as- The sign of Gg cannot be determined in a model-independent way; 
however, it can be predicted by the partonic Hamiltonian in employing naive factor- 
ization. By doing so we find sgn(G8) < confirming the positive interference between 
direct- and indirect- CP- violating components. 



Total rate and present bounds on ^Xt 

Employing the positive sign in (jHHj) we finally arrive at the prediction 

B{Kl ^ 7r°e+e-)sM = {5.2tli) x lO-^^ (42) 

where the error is completely dominated by the uncertainty in B{Ks — > 7r°e'^e~) and, in 
view of (j23); we have neglected the CP-conserving term. The possible improvements of 
this prediction with better data on the Ks mode, and the related sensitivity to '^Xt, are 
illustrated in Fig. IHl 

At present, Eq. ()42j) has to be compared with the new preliminary upper limit by 
KTeV 

B{Kl 7r°e+e-) < 2.8 x 10"^° at 90%C.L., (43) 

based on the full 1997+1999 data sample [6jj. Although we are still far from the possibility 
of precision SM tests in this mode, the combination of the preliminary results by KTeV 
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Figure 3: SM Prediction for B{Ki — > vr'^e+e^) as a function of QXt, neglecting CPC con- 
tributions and assuming a positive interference between direct- and indirect-CP-violating 
components (see text). The three curves correspond to a central value as = 1.08 and 
no error (central full line); 5% error (dashed blue lines); 10% error (dashed green lines); 
present error (red dotted lines). 

and NA48/1 allows us to derive significant bounds on realistic non-standard scenarios (see 
e.g. Ref. |lOj). Instead of discussing any specific model in detail, we can simply express 
these bounds as the limits on dictated by — > 7r°e+e~. In particular, we find 

- 1.3 X 10"^ <Q\t<1.0x IQ-'^ at 90 %C.L., (44) 

which would reduce to \QXt\ < 1-3 x 10^'^ in the absence of any assumption about the 
interference term in (jHEI)- Together with the complementary and comparable constraints 
on Xt derived by — >■ ir'^uu jSSj, these limits represent the most precise information 
presently available about the short-distance structure of the FCNC s —>■ d amplitude. 



5 Dalitz-plot analysis 

The most serious problem in the extraction of the 7r°e"'"e~ amplitude from a 

time-independent rate measurement is the large irreducible background generated by the 
process — >■ 77e+e~ Imposing the cut \M^^ — m^o\ < 5 MeV on the two-photon 
invariant mass spectrum of — > 77e^e~, the latter turns out to have a branching ratio 
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~ 3 X 10 , more than 10^ t imes larger than the signal. As discussed by Greenlee j^I, 
employing additional cuts on various kinematical variables it is possible to reduce this 
background down to the 10^^° level [3 E], but it is hard to reduce it below this figure 
without drastic reductions of the signal efficiency. We stress, however, that this does not 
imply that the signal is unmeasurable in a high- statistic experiment, where the physical 
background can be measured and modelled with high accuracy [35j. 

An important point to notice is that the kinematical analysis necessary to suppress 
and control the Greenlee background provides also a powerful tool to discriminate the CP- 
conserving component of the Kl — > 7i^e~^e~ rate. As discussed in Section El theoretical 
arguments suggest that the CP-conserving component of Kl TT^e~^e~ is very small; 
however, it is clearly desirable to cross-check this statement a posteriori, in a model- 
independent way, using — >■ 7r°e"'"e~ data. 

As discussed by Greenlee the most convenient kinematical variables to describe 
the decay distribution of Kl Ti^e^e~ are 

whose uncorrelated boundaries (in the limit me/ nix 0) are given by 

< 2 < (1 - r^f -l<y <l . (46) 

In terms of these variables, CP-conserving and CP-violating distributions assume the 
following simple factorized structure: 



;cpv 



dz dy 



oc 



y')\'l\l,z,rl)\W{z)\' , (47) 



^'^^""^ - y\l-y')>^ni.z,rl)\B{z)G{z)\^ . (48) 

In Fig. we plot the two distributions, as obtained by setting W{z) = 1 + zm\/my and 
approximating B{z)G{z) to a constant: the y dependence clearly provides a powerful tool 
to discriminate the two terms. 

To better quantify to which extent the two distributions can be distinguished, in Fig. |S1 
we show the result of a cut on y {\y\ < |y|max) on the integrated rate. The three curves 
correspond to the three non-interfering contributions of CP-conserving and CP-violating 
amplitudes, and Greenlee background (as obtained with a constant Kl 7*7* form 
factor), all normalized to one. As can be noted, setting the cut |^|max = 0.5 the CP- 
violating rate is reduced only by 30%, while the CP-conserving one drops almost by a 
factor of 4. The cut on |y|max in not particularly efficient against the Greenlee background, 
but also in this case it improves the signal/background ratio. 



6 Time-dependent interferences 

Complementary information on the direct CP-violating component of the Kl — > iT^e'^e~ 
amplitude can be obtained by studying the time evolution of the Kl^s 7r°e+e~ decay 
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Figure 5: CP- violating (full line) and CP-conserving (dashed line) contributions to the 
Kl —>■ Tc^e^e^ rate (with arbitrary normalization) as a function of |y|max- The dot-dashed 
line denotes Greenlee background. 
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Figure 6: Time dependent distributions for lm{V*gVtd) = (green full), +1.3 x 10""^ (red 
dotted) and —1.3 x 10"'* (blue dashed). Left: probability distribution for a \K'^) state 
at t = to decay into |7r°(77)e"'"e~), Greenlee background included. Right: Probability 
distribution for $ — > K^, Ks — > tt+tt^, 7r°e"'"e~ events at a $ factory (see text). 



^EH]- Although challenging from the experimental point of view, this method has two 
intrinsic advantages: i) the interference between Ks and Kl amplitudes is only due to 
the CP-violating part of the latter; ii) the process Ks — > 77e'^e~ is very suppressed with 
respect to Ks — > vr^e'^e" [B{Ks — > 776+6") ~ few x 10~^^], so the background due to 
the |776"'"6~) final state is almost negligible at small times {t -C tl)- As representative 
examples of time-dependent observables, we shall discuss in more detail two specific cases: 
i) the time evolution of an initial \K^) state, representative of a possible fixed-target 
experimental set up; ii) the time evolution of the KsK^ coherent state produced at a $ 
factory jHTj. 

In the case of a pure beam at t = 0, with Nq particles, the number of decays into 



the final state |7r 6"'"6 ) or in the background channel 
function of the proper time t, can be written as 



"77) (with M^^ ~ m^o) 



as a 



dN 



No 



dt 



m = Ts— = Ts^ \As\' e-'l^' + 2Re A'T A 



-t(rL+Ts)/2TLTs 



+ 



A 



CPV 



+ 



A 



CPC 



+ 



A 



bkg. 



-t/TL 



(49) 



where the decay amplitudes As and Ai^^ are those in Eqs. (j2HI) and PT|) . and an overall 
phase-space integral is understood. The three curves in Fig. IHl have been obtained for 
A^o = 1, assuming ts \As\^ = B{Ks 7r06+6-) = 6 x lO'^ and n l^f^f = B{Kl 
776~^6~)cuts = 10~^°, and employing the following three values of Im(V^*Vj(i): 0, ibl.3 x 
10""^. As can be seen, the interference term is quite sensitive to the value of the direct 
CP-violating amplitude. On a purely statistical level, in this example one could reach a 
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~ 15% error on lm(y^*Vtd) with an initial flux of lO^^/e^Oge -^°5 where e^Oee denotes the 
efficiency for decays occurring within the first 15 Ks decay lengths. 

At $ factories we can take advantage of the quantum properties of the KsKl state 
produced by the $ decay: 

A[(t>-- h{ti)j2{h)] = ^[As{hM)AL{f2M)-As{f2M)AL{hM)] , (50) 

where and ti denote respectively final states and decay times. Choosing /i = Itt+tt^), 
/2 = |7r°e"'"e~) and integrating over ti + t2 [37j, we obtain for t = ti — t2 

(1 + Ts/Tl) 

2t{rL-Ts)/2rLTs , ^5 B{Kl ^ 7r°e + e~) .(^^^_^^^.)/2r,,r.. 



Re 



/H — - 5" L 



(51) 



Also in this case we plot in Fig. [Ul three curves corresponding to Im(V^* V^^) = 0, and 
±1.3 X 10-^ obtained for ts\As\'^ = B{Ks n^e+e-) = 6 x 10"^ and No = 1 but, 
for simplicity, here we ignore the Greenlee background. The direct-CP-violating effect 
appears to be remarkably clear, but the statistics necessary to measure it is probably 
beyond the near-future possibility of existing facilities. 



7 Conclusions 

Motivated by recent experimental results from NA48 on Kl — » ir^'yy and Ks Ti^e^e~ , 
which has been observed for the first time, we have presented a re-analysis of the decay 
mode Kl 7i^e~^e~ . 

The results on Kl —>■ Tc^'~f'~f together with a more general treatment of the Kl —>■ 7r°7*7* 
form factor have led to an updated limit on the CP conserving contribution: B[Kl — 
7r°e"'"e~)cpc < 3 x 10~^^, which we consider as a conservative upper bound. This estimate 
implies that the CP conserving part is essentially negligible for B{Kl 7r^e^e~). 

The new measurement of B{Ks vr^e+e^) has enabled us to fix the size of the indirect 
CP-violating component in Kl — > 7r°e^e~. We have also provided arguments that strongly 
indicate a positive interference with the direct CP-violating amplitude. Combining both 
contributions we predict 

B{Kl ^ 7r°e+e-) = (3.2+^;^) x 10"^^ (52) 

in the Standard Model. The largest contribution is indirect CP-violation, but direct CP- 
violation is also substantial and amounts to ~ 40% of the rate, mostly from interference 
of the two components. The sizable uncertainty which affects at present this prediction 
reflects the poor experimental knowledge of the Ks n^e~^e~ transition. With more 
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precise data on the latter, the theoretical error on B{Kl — >■ ir^e'^e^) could in principle 
be reduced below the 10% level. In this perspective, we stress that a model-independent 
confirmation of the positive interference between direct and indirect CP-violating ampli- 
tudes, which could in principle be obtained by means of lattice QCD, would also be very 
useful. 

As we have briefly discussed, the total rate is not the only interesting observable in this 
channel: the Dalitz plot analysis and the K^-Ks interference, in time-dependent distri- 
butions, could both be very useful in order to extract the direct CP-violating component 
of the Kl — > 7r°e"'"e~ amplitude. 

In summary, our analysis demonstrates that the rare decay Kl — >■ 7r°e"'"e~, despite its 
complexity, can be well described exploiting additional experimental input. It implies, in 
particular, that its branching fraction is larger than might have been anticipated and has 
a substantial sensitivity to direct CP violation and New Physics in AS* = 1 transitions. 
If the experimental challenges posed by a measurement of B{Kl 7r°e+e^) can be 
overcome, this decay will provide us with most valuable information on quark flavour 
physics. 
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Appendix: 0{p^) local contributions to (7^*^)77 

Recent precise data by NA48 [3^ and KTeV |38J, and some recent theoretical papers 
j39| l4Uj ■ have raised an interesting discussion about the determination of the 0{p^) local 
contributions to Kl —>■ 71^77 and the related role of vector-meson contributions [H] . As we 
shall show in the following, a new important piece of information in this respect is provided 
by the very precise measurement of the Kg — ^ 77 rate |12|- Using chiral symmetry, we 
can relate counterterm contributions to Kl — ^ 7r°77 and Ks — > 77 amplitudes. This leads 
to a substantial reduction of the allowed parameter space for the 0{p^) couplings and, as 
a consequence, to a better understanding of the vector-meson-dominance (VMD) ansatz. 

The (8l,1r) weak chiral Lagrangian of 0{p^) contains a huge number of operators. 
However, as long as we are interested only in Kl —>■ tc^'J'J and Kg —>■ 77 decays, we 
can restrict the attention only to three independent combinations ^Hl- In particular, 
we can parameterize all the contributions in terms of the following simplifled effective 
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Lagrangian^ 



47r 



(A.l) 



where, following the notation of Ref. [SI], we define A = uX^u^ 



u, 



iu^D^Uu^ and 



u 



^Mv) + uMu, with U 



u 



e*^''^/^'' {(f) denotes the pseudoscalar-meson field). Ft, 



93 MeV and M. = diag(m^, m^, 2mj^ — m^). The first operator in ()A.1|) contributes both 
to Kl — > 7r°77 and i^'^ 77, whereas the other two can contribute only to Ki —>■ n^jj. 

Historically, the 0{p^) local contributions to A and B amplitudes of Kl — > 7r°77 are 
parameterized in terms of the three unknown coefficients ai, 02 and /3, defined by [T^ 



Act = 0!i{z - rl) + a2, B. 



CT 



/3. 



(A.2) 



Using the Lagrangian ()A.1|) we find 



02 



is: 



i?2 



(4a2 + 203) 



[8ai — 4a2 + 203 



-4a,0 -0.13 



0.65 



(A.3) 



where the extra numerical pieces in 02 and (3 are residual polynomial parts of the loop 
amplitudes, renormalized in the minimal subtraction scheme at the scale = ntp ^HJ- It 
is then easy to check that the VMD ansatz of Ref. ^H] corresponds to the following choice 
for the a,-: 



Oi = 



F3 



■as 



F3 



a2 = 2a 



V- 



(A.4) 



The local contributions in ()A.2|) have to be added to the unitarized loop amplitudes 
to obtain the full result. In the case of the B amplitude this leads to [1(^ 



Biz)=P + cx\^FU- 



2 z' 
+ IO--7 




2, /mr 



(A.5) 



where c is the coefficient, determined from K ^ "in quadratic slopes, which rules the 
strength of unitarity corrections (in the numerical analysis we shall employ the value 
c = 1.1) and the explicit expression of F{z) and R{z) can be found in fTBI. Taking into 
account that F{z) -z/l2 and R{z) z/60 for 2 ^ 0, it follows that 5(0) = -1.71+/3. 
As shown in Eqs. ()14|) - (fT^ . the magnitude of -B(O) is determined in a model-independent 



^ The strongest simplification arises from the fact that we are dealing only with neutral fields, which 
commute with the charge matrix Q. For this reason, the latter is not explicitly inserted in the effective 



operators in l|A.l(l . 
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way by the measurement of B{Kl — ^ 7r°77) at low diphoton invariant mass. Starting 
from the NA48 result 

B{Kl ^ 7r077)M^^e[3o^no] mcV, \ymo-0.2] < 0.6 x IQ-^ at 90%C.L., (A.6) 

which under the assumption of a constant B amplitude becomes 

-> 7r077)M,,<iio Mev < 0.9 X 10-^ at 90%C.L., (A.7) 

we can derive a first clear bound on the 0{p^) counterterms, namely 

2 

-0.4</?<3.8, -1.0 < < 0.07 . (A.8) 

A more precise constraint is obtained by means of B{Ks 77). In this case most 
of the unitarity corrections are implicitly taken into account by the overall coupling Gg, 
extracted from the measured Ks — 7r~^7r~ rate. Non-trivial 0{p^) terms are expected 
only from the (non-VMD) local contributions proportional to Oi, and from the tiny uni- 
tarity corrections associated to the process Kg — > ti^tt^ — 77 [201 • The recent precise 
measurement of NA48, B{Ks 77) = (2.78 ± 0.072) x IQ-^ is substantially higher 
than the 0{p^) prediction, B{Ks — > 77)*-^'' = 2.1 x 10~^ [121, indicating for the first time 
the need for these 0{p^) terms. Summing the local term generated by Cq to the leading 
loop amplitude of 0{p'^) j321 we can write 



A{Ks ^ 77) = T^iMl - M] 



irm 




F2 «^ 



^i^^2u [qUi - {<li<l2)9 



(A.9) 

where the relative sign between local and non-local terms has been unambiguously fixed by 
the sign convention adopted in ()A.3|1 . Using oi as a free parameter to fit the experimental 
branching ratio we finally obtain 



5m^ 



2 

ai = {ai + a2 + (3) + 0.78 = 1.0 ± 0.3 . (A.IO) 



The error in ()A.in|l is not due to the experimental uncertainty, but is a theoretical estimate 
of possible subleading terms: the ±lcr value in ()A.lfl|l corresponds to the following two 
cases: i) Eq. flA.9|) only; ii) Eq. ()A.9|) plus the absorptive contribution due to Ks 

7r°7r° —>■ 77. 

Taking into account the constraints in ()A.10|) and ()A.8|1 we are basically left with 
a single free parameter to fit both rate and high diphoton invariant mass spectrum of 
Kl 7r°77. Apparently, the constraint on jS in ()A.8|) is not very stringent; however, 
this does not represent a serious limitation. The situation becomes much more clear if we 
adopt the S-B basis for the Kl — ^ 7r°77 amplitude {S = A + B): in this framework the 
bound in ()A.8|) tells us that the B amplitude plays a very marginal role - but for very low 
diphoton invariant masses (see Section|2I) - and that both rate and high diphoton invariant 
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mass spectrum of Kl 7r°77 are completely dominated by S. The latter depend only 
on two independent counterterm combinations, namely ai and (02 + P), whose sum is 
severely constrained by ()A.10|) . It is then clear that we are left with a single effective 
coupling, which we chose to be ai. Summing the local contributions to the unitarized 
loop amplitude and fitting the NA48 result B{Kl vr°77) = (1.36 ± 0.05) x 10"*^ , we 
find 

2 

«i = ^(4a2 + 2a3) =3.4±0.4 , (A.ll) 

where the error reflects mainly the uncertainty in ()A.10|) .^ 

A detailed discussion of the results in ()A.8|) . ()A.10|) and ()A.11|) goes beyond the scope 
of this work and will be presented elsewhere. The quoted figures should also be taken 
with some care, since they are not based on a complete analysis. On the other hand, 
we can already draw a few interesting conclusions: 

• The local structure of Ks — > 77 and Ki — > tt" (77)7=9 amplitudes is strongly con- 
strained by chiral symmetry. Taking into account this theoretical constraint and the 
recent precise data by NA48, we are able to fix in a precise way the two indepen- 
dent combinations of counterterms, ai and 202 + 03, which control these amplitudes. 
Interestingly, these results rule out completely or strongly constrain most of the so- 
lutions proposed in [IH] , taking into account Kl — > 71^77 data only. 

• As discussed in Section 01 the smallness of the Kl (77)7=2 amplitude is a 
direct consequence of the experimental result in ()A.6|) . We thus agree with the 
conclusion of Ref. jlUI that the estimate of the CPC contribution to Kl —>■ 7r°e+e~ 
does not depend on specific assumptions about Kl — ^ 71^77 counterterms. Since 
the low diphoton invariant mass analysis of NA48 contradicts previous findings by 
KTeV [38j, an independent confirmation of ()A.6|) would be very welcome. 

• The evidence for a non-vanishing Oi in ()A.10|) shows that the VMD ansatz ()A.4|1 
is not exactly fulfilled. Nonetheless, VMD contributions are likely to provide the 
dominant 0{p^) contribution. This statement is not very quantitative at the mo- 
ment due to the uncertainty on 03 in ()A.8|) . However, we note that the values for 03 
preferred by the low diphoton invariant mass analysis of NA48 are those close to the 
boundaries of ()A.8|) : {m\ / F'^)az ^ —0.9 or 0.1. In other words, the B amplitude is 
likely to be just below the exclusion limit implied by ()A.6|1 . If (m|-/F^)a3 ~ —0.9, 

® The central value in (|A.lip has been obtained with the inclusion of the absorptive contribution due 
to Kl — > Stt" 7r''77 in the non-local A amplitude jSH): with this choice, we find that the predicted 
Kl 7r°77 spectrum is in excellent agreement with the observation of NA48 [Sj- If the absorptive 
contribution due to Kl —>■ Stt" —> tt'^jj is not included, the central value of ai fitted from the rate turns 
out to be about 20% higher, but the resulting diphoton spectrum is not in very good agreement with 
NA48 results. We have also explicitly checked that the central value in (|A.11I) is completely insensitive 
to the choice of f3, varied within the interval l|A.8|l . 



21 



taking into account also the constraint ()A.11|) . we would conclude that 




{0,3 - 02) 



(as + ^2) 



+0.4 



2.2 



VMD 



VMD 



4av , 



(A.13) 



(A.12) 



which can be interpreted as a clear manifestation of vector dominance. 
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